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ABSTRACT 

We report on HST observations of six candidate old globular clusters in the Large Mag- 
ellanic Cloud: NGC 1754, NGC 1835, NGC 1898, NGC 1916, NGC 2005 and NGC 
2019. Deep exposures with the F555W and F814W filters provide us with colour- 
magnitude diagrams that reach to an apparent magnitude in V of ~25, well below 
the main sequence turnoff. These particular clusters are involved with significantly 
high LMC field star densities and care was taken to subtract the field stars from the 
cluster colour-magnitude diagrams accurately. In two cases there is significant variable 
reddening across at least part of the image, but only for NGC 1916 does the differ- 
ential reddening preclude accurate measurements of the CMD characteristics. The 
morphologies of the colour-magnitude diagrams match well those of Galactic globular 
clusters of similar metallicity. All six have well-developed horizontal branches, while 
four clearly have stars on both sides of the RR Lyrae gap. The abundances obtained 
from measurements of the height of the red giant branch above the level of the hori- 
zontal branch are 0.3 dex higher, on average, than previously measured spectroscopic 
abundances. Detailed comparisons with Galactic globular cluster fiducials show that 
all six clusters are old objects, very similar in age to classical Galactic globulars such 
as M5, with little age spread among the clusters. This result is consistent with ages 
derived by measuring the magnitude difference between the horizontal branch and 
main sequence turnoff. We also find a similar chronology by comparing the horizontal 
branch morphologies and abundances with the horizontal branch evolutionary tracks 
of Lee, Demarque, & Zinn (1994). Our results imply that the LMC formed at the same 
time as the Milky Way Galaxy. 

Key words: Magellanic Clouds - stars: Population II - galaxies: star clusters - 
galaxies: formation. 



1 INTRODUCTION 

It has long been known that the LMC contains several clus- 
ters that strongly resemble normal globular clusters of the 
Galaxy (Shapley 1930). Several of these objects, however, 
have turned out to be of intermediate or even young age, 
in spite of their outward appearance of large size and lu- 
minosity. Some years ago one of us, in an attempt to use 
crude colour-magnitude diagrams, selected 35 clusters that 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc. under NASA contract NAS5-26555. 



seemed to be old globular clusters (Hodge 1960). Truly old 
clusters are actually rarer, however; a recent review of the 
old populations of the MCs lists only 14 objects that are 
candidates for clusters that might be as old as the oldest 
Galactic globulars (Olszewski, Suntzeff and Mateo 1996; see 
also Westerlund 1997). However, precise ages are very dif- 
ficult to determine from the ground, as the main sequence 
turnoff (hereafter MSTO) occurs at a y of approximately 23, 
too faint for accurate ground-based photometry in crowded 
regions of the LMC. For the less crowded outlying clusters 
NGC 1466, NGC 1786, NGC 1841, NGC 2210, NGC 2257, 
and Reticulum, ground-based CMDs approaching and going 
below the MSTO have demonstrated that they are quite old 
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Table 1. Basic LMC Globular Cluster Parameters 



Cluster 


RA (2000.0) 


Dec 


V (int.) 


n(RR) 


NGC 1754 


04h55m 


-70?31' 


11.4 




NGC 1835 


051i05m 


-69? 28' 


9.5 


35 


NGC 1898 


05hl7m 


-69? 43' 


11.1 




NGC 1916 


05hl9m 


-69? 27' 


9.9 




NGC 2005 


05h30m 


-69? 45' 


11.2 




NGC 2019 


051i32m 


-70? 12' 


10.7 






(Brocato et al. 1996; Walker 1990, 1992a, 1992b), with NGC 
2257 perhaps being ^2 Gyr younger than the bulk of Milky 
Way clusters (Testa et al. 1995) . With uncertainties in age of 
typically ^3 Gyr, however, the data are not accurate enough 
to fully determine whether the clusters are truly as old as the 
classical Galactic globulars or are 2-3 gigayears younger, as 
has been suggested on the basis of their horizontal branch 
(hereafter HB) morphologies (Da Costa 1993). This ques- 
tion is of strong interest both for establishing the timescale 
of formation of the LMC and for testing age as the candi- 
date second parameter affecting IfB morphology, a topic of 
currently hot debate (cf. Lee 1992; Lee, Demarque, & Zinn 
1994, hereafter LDZ; Stetson, VandenBerg, & Bolte 1996; 
Sarajedini, Chaboyer, & Demarque 1997). 

This paper reports on a Hubble Space Telescope pro- 
gram to determine colour-magnitude diagrams for the six 
inner globular clusters of the LMC, for which ground-based 
data would be especially difficult to interpret because of the 
photometric effects of crowding. The clusters and their char- 
acteristics (taken from Olszewski et al. 1996) are listed in Ta- 
ble 1. These clusters, which presumably formed deep in the 
LMC gravitational potential well, are excellent probes of the 
early formation of the main mass of the LMC. In conjunction 
with the HST study of Hodge 11 by Mighell et al. (1996) and 
the comprehensive study with HST of the outlying old LMC 
clusters currently underway (Johnson & Bolte 1997), our re- 
sults will be useful for testing models of the formation of the 
LMC, much as observations of globular clusters in the Milky 
Way have evaluated formation models of the Galaxy. In the 
following we describe the observations in Section 2, detail 
our reductions in Section 3, discuss field star decontamina- 
tion in Section 4, and derive the colour-magnitude diagrams 
(CMDs) for the clusters in Section 5. In Section 6, we present 
the ages of the clusters implied by three age-dating tech- 
niques: comparison of our CMDs with Milky Way globular 
cluster fiducial sequences through the "horizontal method" 
(VandenBerg, Bolte, & Stetson 1990, hereafter VBS), mea- 
surement of AV^HB ("vertical method"), and comparison of 
our HB data with HB evolutionary models from LDZ. These 
ages are dependent on the adopted abundances, which we 
measure from the CMDs and compare with previously mea- 
sured spectroscopic abundances (Olszewski et al. 1991, her- 
after 091). We also describe in Section 6 our estimates of 
the reddenings and distances of the clusters. These estimates 
are not direct measurements, but rely on knowledge of the 
reddenings and distances to Galactic globular clusters and 
on an assumed Mv(RR)-[Fe/H] relationship. 



2 OBSERVATIONS 

Observations of the clusters were taken during Cycle 5 of 
HST. Each cluster was observed for the duration of one orbit 
of the spacecraft. Table 2 contains a summary of the observa- 
tion log. For these observations, we used the low gain setting 
of the WFPC2 camera (GAIN=7). Multiple exposures were 
taken through each filter to aid in cosmic ray removal, but no 
dithering technique was used. We took both long and short 
exposures to provide unsaturated photometry of as many 
stars as possible. The long exposures reach limiting mag- 
nitudes of ~24.5 in F555W and ~23.5 in F814W, or ~1.5 
magnitudes below the main sequence turnoff in V . The im- 
ages were processed through the standard STScI reduction 
pipeline prior to our receipt of them. 

In order to check the zero point of the HST photometry, 
one of us (A.W.) obtained images of each cluster with the 
CTIO 1.5-m and Tek 2048 CCD. The field of view of the 
camera is 8.'3, or ~3 times the size of the WFPC2 field 
of view. Images were taken through CTIO copies of the 
//SrF555W and F814W filters on the nights of January 23- 
26, 1995 under photometric conditions and in good seeing 
(FWHM~1"). Standard stars from Landolt (1992) and from 
the u) Gen field used to calibrate WFPC2 photometry were 
observed. Careful photometry using DAOPHOT/ALLSTAR 
was performed on suitably isolated stars in the ground-based 
frames and the photometry transformed to Johnson/Kron- 
Cousins V and I. The comparison of the photometry from 
these ground-based images with the HST photometry is de- 
scribed in section 3.6. 



3 REDUCTIONS, PHOTOMETRY, AND 
CALIBRATION OF WFPC2 IMAGES 

3.1 Image reduction and corrections 

Our first step was to check the image alignment for each 
set of exposures in each filter. In all cases, the alignment is 
better than a fraction of pixel. We made no attempt to align 
the images further. 

We removed the cosmic rays from the images with the 
STSDAS task crrej. From multiple equal-length ex- 
posures of the same field in each filter, crrej produces a cos- 
mic ray-free image by averaging the stack of exposures after 
rejecting pixels with values too high compared to an ini- 
tial guess at the uncontanimated pixel values. The threshold 
outside of which pixels are rejected is set by the noise char- 
acteristics of the image and through adjustable parameters. 
Because there are fractional pixel ofi'sets between successive 
frames, the sharp cores of stars illuminate the pixels of the 
successive frames difi'erently, and are often interpreted as 
cosmic rays by the standard rejection procedure. By setting 
the parameter scalenoise to 10 per cent, which raises the 
rejection threshold by 10 per cent of the pixel value, we in- 



t IRAF is written and supported by the IRAF programming 
group at the National Optical Astronomy Observatories (NOAO) 
in Tucson, Arizona. NOAO is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. under cooperative agree- 
ment with the National Science Foundation. 
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Table 2. Summary of Observing Log 



Cluster 


Files 










Filter 


Exp. time (sec.) 


Date (DD/MM/YY) 


i\ V_l V_/ _L ( tj'-t 






ii9xnnin9t rvt 






J- <J'J<J V V 


20 


21 /in/Q'i 




u2xq0103t cvt 




u2xq0104t cvt 




u2xq0105t cvt.* 


F555W 


500 


21 /1 0/95 




ii9-x-nni Dfit ri't 




ii9xnnin7t rvt 




n9yn01 rvt * 


F814W 


20 


21 /in/95 








ii9xnninat rvt 




n9vnninht rvt * 


F814W 


600 


21 /in/95 








ii9-ifnn9n9t rvt 






F555W 


20 


18/10/95 








n9vnn9n4t rvt 




n9vnn9n'St rvt * 


F555W 


500 


18/10/95 








n9-5fnn9n7t rvt 




n9vnn9nSt rvt * 


F814W 


20 


18/10/95 




n9vnn9nQt rvt 




ii9Ynn9nat rvt 




n9vnn9nht rvt * 


F814W 

L O J-TT V V 


600 


18/10/95 


NGC 1898 


ii2xnn801t rvt 




Ii2xo0302t rvt 






F555W 


20 


10/12/95 

x\^/ x^/ c 




ii^vnO'^O'^t rvt 




u2xq0 304t _c vt 




u2xq0305t_cvt. * 


F555W 


500 


10/12/95 

XW/ X^/ T_7^_J 




u2xq0306t_cvt 




u2xq0307t_cvt 




u2xq0308t_cvt. * 


F814W 


20 


10/12/95 

xu/ x^/ c 'J 




u2xq0309t_cvt 




u2xq030at cvt 




u2xq030bt cvt.* 


F814W 


600 


10/12/95 

xu/ x^/ c f-j 


NGC 1916 


u2xq0401t_cvt 




u2xq0402t_cvt 






F555W 


20 


10/12/95 

xu/ x^/ c o 




ii2xn0403t rvt 




ii9xnn4n4t rvt 




ii9ynn4.0^t rvt * 


F555W 


500 


10/12/95 




n9vnn4nfit rvt 




n9-5fnn4n7t rvt 




n9vnn4n8t rvt * 


F814W 


20 


10/12/95 




u2xq0409t.cvt 




u2xq040at_cvt 




u2xq040bt.cvt.* 


F814W 


600 


10/12/95 


NGC 2005 


u2xq0501t.cvt 




u2xq0502t_cvt 


* 




F555W 


20 


19/10/95 




u2xq0503t.cvt 




u2xq0504t_cvt 


*, 


u2xq0505t_cvt.* 


F555W 


500 


19/10/95 




u2xq0506t_cvt 




u2xq0507t_cvt 


*, 


u2xq0508t.cvt.* 


F814W 


20 


19/10/95 




u2xq0509t_cvt 




u2xq050at_cvt 


*) 


u2xq050bt_cvt.* 


F814W 


600 


19/10/95 


NGC 2019 


u2xq0601t_cvt 




u2xq0602t_cvt 


* 




F555W 


20 


18/10/95 




u2xq0603t_cvt 




u2xq0604t_cvt 




u2xq0605t_cvt. * 


F555W 


500 


18/10/95 




u2xq0606t.cvt 


*) 


u2xq0607t_cvt 




u2xq0608t.cvt.* 


F814W 


20 


18/10/95 




u2xq0609t_cvt 


*5 


u2xq060at_cvt 




u2xq060bt_cvt.* 


F814W 


600 


18/10/95 




Figure 1. Fit of the DoPHOT analytical PSF used in the photometry to a sample star from the WFPC2 frames. While the PSF fits 
the core of the star and the background well, it doesn not fit the broad wings most prominent in the F814W images. 



sured that the cores of stars were left unaltered while still 

removing many cosmic rays. 

Wc accounted for CTE effects by using the corrections 
suggested by Holtzman et al. (1995a). The backgrounds in 
the long exposures are in the 30-200 e~ range, while the 
short exposures have backgrounds of 2-10 e~'s. For the short 
exposures, wc multiplied an image containing a 4 per cent 
ramp into the data frames, while we used a 2 per cent ramp 
for the long exposures. Subsequent to our image reduction, 
Whitmore & Heyer (1997) published a report indicating that 



CTE corrections should be a function of x- and y-position 
on the chip as well as star counts and background counts, 
but the effects are small and we have not included them. 

Because we are doing point source photometry, we need 
to correct for the fact that the flat-fielding process adjusts 
the pixel counts so that surface brightness is preserved with 
position on the chip. We used the Holtzman ot al. (1995a) 
coordinate transformations to create pixel distortion maps 
for each chip, which we multiplied into the data frames. 

We used the DQF images to create bad pixel masks for 
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each group of exposures in each filter. These maslis include 
macroscopic charge traps, bad columns, questionable pixels, 
and pixels saturated by stars in the field. By sotting the 
masked pixels in the data frames to a low value, they are 
ignored during photometry. 

3.2 DoPHOT Photometry 

In order to achieve the best photometry for the faint stars, 
which are crucial for determining cluster ages, it is necessary 
to resort to profile-fitting rather than aperture photometry. 
For its speed and ease of use, we chose to use the program 
DoPHOT (Schcchtcr, Mateo, & Salia 1993), version 2.5, with 
modifications performed by Eric Deutsch to handle our float- 
ing point images. DoPHOT's operation is controlled by a file 
containing user-adjustable parameters, the most important 
of which are described in detail below. 

Wc chose to use DoPHOT's median filter to generate 
the model background above which DoPHOT identifies ob- 
jects. We found that DoPHOT's other background options, 
the plane and Hubble models, did not adequately fit the 
background, becoming weighted too heavily by the bright 
cluster cores and causing stars at the outskirts of the clus- 
ters to escape detection. We used a minimum threshold for 
object detection of 2a above the background. 

We adjusted the coefficients of DoPHOT's power law 
point-spread function (PSF) by simultaneously inspecting 
the fit of the PSF to a small sample of stellar radial pro- 
files and by minimizing the aperture correction, described 
in Section 3.3. We fit the coefficients independently for each 
filter and PC or WF chip. Table 3 contains the final adopted 
values of 134 and /Se, and Fig. 1 shows a sample fit of the PSF 
to a stellar radial profile. Although the PSF fits well in the 
core of the star, the wing shows structure that is not fit by 
the PSF, most noticeably in the F814W filter. Hence, the 
PSF-subtractcd imago shows halos around the positions of 
subtracted stars. In order to prevent DoPHOT from identi- 
fying too many spurious objects in these residuals, we added 
extra noise to DoPHOT's noise array. 

Finding the optimal size of the box within which the 
PSF is fit is a trade-off between the need to include pixels 
far from the core of the star for a proper fit to the back- 
ground and the desire to exclude neighboring stars from the 
fit box. The issue of crowding on the accuracy of the pho- 
tometry is most critical in the crowded PC frames. Fig. 2 
shows the distribution of sky levels with position on the PC 
chip for two fit box sizes, 9 and 25 pixels, and the resulting 
CMDs of NGC 1754. While the use of the larger 25-pixel 
fit box produces a much smoother sky distribution, the ac- 
companying CMD shows much greater scatter. We settled 
on the 9-pixel fit box for both the PC and WF frames, find- 
ing that this size produced adequate determinations of the 
background while avoiding, as much as possible, the effects 
of crowding. 

When DoPHOT encounters objects that are signifi- 
cantly broader than the mean stellar profile, it makes a de- 
cision on whether to interpret the object as a single broad 
object ("galaxy") or as overlapping stars, in which case the 
object is split into multiple objects. The decision is weighted 
by the parameter stargalknob, which we set to a high value 
because these crowded frames likely contain many more 
overlapping stars than background galaxies. We encountered 



Table 3. DoPHOT PSF Parameters 

F555W F814W 

Chip /34 ft FWHM Pi pe FWHM 

PC 4.0 1.5 2.0 4.5 4.5 2.0 

WF2-4 2.2 2.2 1.2 2.5 2.5 2.0 



some problems with DoPHOT classifying some bright yet 

unsaturated stars in the long exposure frames as broad, and 
splitting them into pairs. Using the short exposure photom- 
etry for these cases solved this problem. 

3.3 Aperture corrections 

Because the model PSF is imperfect, the PSF-fitting pho- 
tometry produces systematically different results from aper- 
ture photometry. This systematic difference varies between 
different chips and filters. Moreover, because the WFPC2 
PSF varies with position on the chip, while DoPHOT as- 
sumes that the PSF is uniform with position, we applied an 
aperture correction that depends on position in the frame 
as well as the chip and filter used. 

For these aperture corrections, it was necessary to find 
several bright, isolated stars in order to be able to do aper- 
ture photometry to the 0'.'5 radius on which the Holtzman et 
al. (1995b) calibrations are based. While we found a few such 
isolated stars in the WF frames, there were practically none 
in the PC frames. We settled for using a PSF to subtract all 
stars but a few bright ones and performing aperture photom- 
etry on the subtracted image. Because DoPHOT's PSF is 
uniform in position and leaves halos in the subtracted image, 
wc chose to use PSFs kindly provided by Peter Stetson and 
DAOPHOT/ALLSTAR to produce the subtracted images. 
These PSFs, derived from multiple WFPC2 observations of 
iu Ccn, arc "perfect" PSFs, in that they average over time- 
dependent effects such as focus variations. We found that 
ALLSTAR photometry using the ut Cen PSFs had greater 
scatter than our DoPHOT photometry, although the PSF- 
subtracted images are cosmetically cleaner. 

Choosing 200 bright stars to be left untouched in each 
frame, we used the ALLSTAR photometry to subtract the 
remaining stars. We inspected the profiles of each of the 200 
stars in the subtracted frame, discarding stars that had sig- 
nificant unsubtracted companions or deviant pixels within 
a 0'.'5 radius. With DAOPHOT's PHOTOMETRY routine, 
we performed aperture photometry on each of the 200 stars 
out to 0'.'5. After matching the aperture photometry lists 
with the DoPHOT photometry list, we fit the aperture cor- 
rections as a function of x and y on the chip with an appro- 
priate polynomial. Figs. 3-4 show sample fits of the aperture 
correction surfaces in the PC and WF3 frames for both long 
and short exposures. In all of the WF frames, a 2nd or- 
der polynomial surface fits well. The residuals of the points 
around the mean surfaces are generally Gaussian-distributed 
with dispersions similar to those expected from errors in the 
combined aperture and PSF photometry, implying that we 
were successful in removing the stars neighboring to those 
used in calculating the aperture corrections. However, in the 
PC frames we find no clear dependence of the aperture cor- 
rection on position, which we interpret as the result of the 
severe crowding and the difficulty of doing reliable aperture 
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Figure 2. Plots showing tlie eflect of diflferent sizes of the fitting box used in the DoPHOT photometry. While using a larger fitting box 
produces better background fits, crowding limits the a<;curacy of the photometry. 



photometry in the PC frames. Instead of using the question- 
able PC aperture photometry, we tested using an aperture 
correction based on Stetson's PSFs. The artificial star tests 
used to generate the aperture correction surfaces based on 
Stetson's PSFs arc described in Section 4. To test the suit- 
ability of using an artificially generated aperture correction 
surface in the PC frames. Fig. 5 shows a sample comparison 
between the WF2 aperture correction surfaces generated us- 
ing our own data with artificial surfaces. In most cases the 
agreement is good, but for some of the frames there are dif- 
ferences between the surfaces of ~0.05 mag. However, as we 
are uncertain of the PC aperture corrections generated from 
our own data at a level of >0.05 mag, we chose to use the 
artificially generated surfaces to correct the PC photometry. 



3.4 Offset between short and long exposure 
photometry 

In agreement with the report of Whitmore & Heyer (1997), 
among others, on CTE effects in WFPC2 photometry, we 
find an offset between magnitudes measured in the long 
exposure frames compared to those measured in the short 
exposure frames. Fig. 6 shows a sample comparison. The 
difference increases as the magnitude of the star increases, 
reaching a maximum of ~0.1 magnitudes at F555W = 21.5. 
As the offsets are due to a CTE effect which worsens at 
low background and count levels, we chose to trust the long 
exposure photometry and apply a correction to the short ex- 
posure photometry. Within broad magnitude bins, we cal- 
culated the average difference between the long and short 
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exposure magnitudes and subtracted this difference from 
all of the short exposure magnitudes in the bin. The use 
of broad magnitude bins assumes that the offset changes 
slowly with magnitude, with no erratic behavior hidden by 
the scatter in the plots. We checked this assumption by com- 
paring the distribution of points around the average offset 
within each bin to the equivalent distribution generated by 
the artificial star tests described in Section 4. The appli- 
cation of a Kolmogorov-Smirnov statistic shows that in all 
cases the distributions arc nearly identical, implying that the 
assumption that the offset changes slowly with magnitude 
is well-justified. 



3.5 Final calibrated photometry 

We merged the short and long exposure photometry lists as 
follows. First, we removed stars near saturated and bad pix- 
els from the lists. Stars within five pixels of a saturated pixel 
were assumed to lie in the wings of the saturated star and 
were removed, while stars within 1.5 pixels of a bad pixel 
wore also removed. Next, wo matched the lists according to 
position of the stars. When we plot the distribution of pixel 
offsets between stars found in both the short and long expo- 
sure lists, we find that most of the distribution lies within 
0.6 pixels of the origin, so matches outside this range were 
considered spurious. For all matching pairs, wo kept the pho- 
tometry for the star having the smallest photometric error. 
Sometimes, a star identified as single on the short exposure 
image was split into two stars by the long exposure reduc- 
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Figure 4. An example showing the much higher quality of the aperture corrections of the less crowded WF frames. The solid lines are 
again 2nd order polynomial fits while the dashed lines are Icr deviations. 



tion. In these cases, we compared the photometric error of 
the single star with the error of the split pair member with 
the position closest to the single star, and kept the data for 
the star having the smallest error. In most cases, this pro- 
cedure rejected the photometry of the split pair member in 
favor of the single star and generally eliminated the previ- 
ously mentioned difficulty of having bright stars split into 
two in the long exposures. We also added stars found in the 
short exposure frames but not in the long exposures to the 
final combined list, as these stars were generally saturated 
in the long exposures. 

Next, we matched the merged photometry lists from 
the two different filters, again using a 0.6 pixel matching 
radius. We kept only stars having a DoPHOT object type of 
1 ("perfect" star), 2 ("galaxy"), 3 (member of overlapping 
stellar group), or 7 (faint but measurable object) for the 



final list. Using equation 9 and the zero points from table 10 
of Holtzman et al. (1995b), we transformed our magnitudes 
to Johnson V/Kron-Cousins /, as this is the photometric 
system on which stellar evolutionary models are generally 
tabulated. Figs. 7a-f show the CMDs generated from the 
PC photometry. These CMDs are dominated by the cluster 
stars, but show some contamination due to field stars. Table 
4 contains a sample of our photometric data, the rest of 
which is avail able electronically through the Astron omical 
Data Center ( ittp : //adc . gsf c .nasa. gov/ adc .html). 



3.6 Comparison with ground-based photometry 

In order to check the Holtzman et al. (1995b) zero points, 
we checked our calibrated WFPC2 photometry against our 
V and / CTIO 1.5-m observations of the cluster fields. The 
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Figure 6. A sample comparison of stellar magnitudes measured in the long (500, 600 s) WFPC2 exposures with the short (20 s) 
exposures, showing the CTE eflfect discussed by Whitmore & Heyer (1997). 



Table 4. LMC Globular Cluster Photometry 



Star # 


X(V) 


Y 


V 


"v 


/ 


a- 


V type' 


I type 


NGC 1754 


















1 


386.11 


100.57 


19.3881 


0.0380 


18.3591 


0.0520 


10 


10 


2 


477.04 


104.77 


18.6639 


0.0420 


17.5141 


0.0490 


3 


10 


3 


491.68 


112.58 


19.6949 


0.0440 


18.6781 


0.0490 


1 


10 


4 


597.60 


127.65 


19.4229 


0.0500 


18.3281 


0.0510 


1 


10 


5 


233.08 


129.43 


19.6629 


0.0460 


19.3249 


0.0540 


1 


1 


6 


142.73 


142.05 


19.6841 


0.0390 


19.5021 


0.0560 


10 


10 


7 


753.78 


153.25 


19.3071 


0.0550 


18.1461 


0.0500 


30 


10 


8 


286.41 


184.42 


19.2699 


0.0460 


18.8031 


0.0450 


1 


10 


9 


105.07 


203.56 


19.6689 


0.0540 


19.2081 


0.0520 


3 


10 


10 


644.67 


212.54 


19.5009 


0.0470 


18.4271 


0.0480 


1 


10 



Removed'^ 



"Photometric errors reported by DoPHOT 

''Where the short exposure photometry was used, the DoPHOT object type has been multiplied by 10. 
See text for explanation of object types. 

•^Stars removed by field star cleaning procedure are marked with "y". 



seeing in the 1.5-m images is generally corresponding 
to 22 pixels on the PC and 10 pixels on the WF chips. Using 
approximate WFPC2 x and y positions of the ground-based 
stars, we identified stars from the 1.5-m and WFPC2 lists 
with positions matching within a radius equal to the ra- 
dius of the 1.5-m seeing circle. For almost all of the stars 
appearing isolated on the 1.5-m images, we found several 
stars within the seeing circle on the WFPC2 frames. The 
correct match was assumed to be with the brightest of the 
candidates in the WFPC2 frame. We inspected each match- 
ing star by eye on the WFPC2 frame, and eliminated any 
star with significantly bright companions within a 0'.'7 radius 
from further consideration. The remaining matches were 
used to establish the comparison between the ground-based 
and WFPC2 photometry. 

Because the four WFPC2 chips are essentially treated 
as different instruments in our reduction, we expect that 



the average differences between the ground-based and HST 
magnitudes, measured on each chip and through both filters, 
will be randomly scattered, independent of chip or filter. In- 
stead, we find that they are strongly correlated. We attribute 
this correlation to the extreme difficulty of doing photom- 
etry in the severely crowded ground-based frames. We also 
find that the difference between the ground-based and HST 
magnitudes is correlated with distance from the cluster cen- 
ter, which we believe is due to the increased crowding to- 
wards the cluster core. After removing the cluster-to-cluster 
variations, the ground-based and WFPC2 zero points are 
found to agree to within ~0.1 mag. We have thus not felt it 
necessary to adjust the Holtzman (1995b) zero points. 



4 ARTIFICIAL STAR TESTS 
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Figure 7. Color-magnitude diagrams of the six old LMC globular clusters studied here. No field star subtraction or other selective 
removal of stars has been performed and the diagrams are uncorrected for reddening. The error bars are derived from the artificial star 
tests discussed in the text. 



4.1 Generating the artificial star lists 

We conducted artificial star tests to study the complete- 
ness of tfie pliotometry, to examine tfie effects of crowd- 
ing, and to aid in separating tlie field stars from the cluster 
stars. The crowding and completeness effects are expected 
to be a function of magnitude, colour, and position in the 
frame. Therefore, wc used a set of artificial stars with mag- 
nitudes, colours, and positions distributed similarly to those 
of the stars in our CMDs. We generated the magnitudes 
and colours of the input artificial star sets by dividing our 



combined WF and PC CMDs into bins of 0.03-0.06 magni- 
tudes inV — I and 0.1-0.2 magnitudes in V. In each bin, we 

chose a number of random points to represent the artificial 
stars. The number of artificial stars in a bin was chosen to 
be 10 times the rmmbor of observed stars in that bin, with 
the constraint that no bin should contain more than 100 
artificial stars. Fig. 8a shows a sample input artificial star 
CMD. We distributed the positions of the artificial stars 
randomly, weighting the distribution by the stellar density 
profile. Wc determined the stellar density profile by fitting 
the completeness-corrected data to a King model plus a con- 
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Figure 8. An example of the input artificial star distribution (a) used in the tests discussed in the text and the colour-magnitude 
diagram recovered after addition to and processing through the WFPC2 images (6). 




Figure 5. An example of the agreement between aperture cor- 
rection surfaces generated using Stetson's PSFs (dashed lines) 
with those derived from our WF2 frames (solid lines). Both the 
2nd order polynomial fits to the surfaces and the la deviations 
are shown. The complete set of comparisons is available in Olsen 
(1998). 



updated completeness corrections of the new artificial star 
experiments. 

Using the previously mentioned PSFs obtained from 
Stetson, we distributed the artificial star sets over multiple 
copies of the original images, with the density of artificial 
stars in a given image amounting to ~5 per cent of the orig- 
inal stellar density. Each artificial star was placed in both 
the long and short exposure frames and in both filters. In all, 
~6000 artificial star images were generated. The frames con- 
taining the artificial stars were reduced in exactly the same 
manner as the original images, although it was, of course, 
unnecessary to perform any CTE correction to the short ex- 
posure photometry. To recover the detected artificial stars, 
we searched the DoPHOT output lists for stars with posi- 
tions matching within 0.6 pixels to those of stars in the input 
artificial star list. The list of recovered stars was similarly 
compared with the list of stars from the original image. In 
cases where the recovered artificial star also matched a star 
in the original image, we considered the artificial star lost if 
the recovered magnitude was closer to the magnitude of the 
real star than to the magnitude of the artificial star. As with 
the real photometry, artificial stars placed too closely to bad 
or saturated pixels were considered lost. We merged the lists 
of recovered artificial stars placed in the short and long ex- 
posure images, and matched the merged lists from the two 
filters. Fig. 8b shows the CMD of the stars recovered from 
the distribution in Fig. 8a. 



stant background, with completeness corrections measured 
from am earlier artificial star experiment in which the arti- 
ficial stars were uniformly distributed in position. Although 
useful for estimating the completeness, this initial experi- 
ment was not helpful in statistically subtracting the field 
and cluster CMDs from one another because of insufficient 
coverage in V — I. We refit the King models later from the 



4.2 The effects of crowding and incompleteness 

Crowding and incompleteness aifect our analysis in a num- 
ber of ways. Because the differences in crowding between the 
cluster-dominated PC frames and the field star-dominated 
WF frames are large, we need to quantify the crowding ef- 
fects in order to separate the field and cluster stars. Proper 
field and cluster star separation also depends on our fit to 
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NGC 1754 




Figure 9. A sample of the calculated completeness as a function 
of V and radius from the cluster center. The gray shadings arc 
levels of constant completeness fraction. At bright V magnitudes, 
the completeness has been set to 1 where no artificial stars were 
placed to calculate the completeness. 



the cluster profile. To measure this profile from the stellar 

density distribution, wc need to measure the completeness 
as a function of both position and magnitude. Finally, be- 
cause the cluster ages that we measure depend sensitively 
on the accuracy of the photometric scale, we need to know 
whether crowding introduces any changes in the colour or 
magnitude scales. Based on our understanding of the crowd- 
ing and completeness through the artificial star tests, we will 
be able to select those stars with the most reliable photom- 
etry for conducting our analysis of the CMDs. 

Fig. 9 shows a sample completeness surface calculated 
as a function of position and V magnitude. The surface rep- 
resents the probability that a star of a given position and 
actual V magnitude would be found in our CMDs. The sur- 
face shows how the completeness and associated limiting 
magnitude are strong functions of position near the cluster 
center. 

We used the artificial star tests to evaluate both random 
errors and systematic colour and magnitude shifts in the C- 
M plane. Fig. 10 shows an example of the random errors and 
systematic shifts in the photometry. While the size of the 
random errors depend mostly on V with very little, if any, 
colour dependence, the systematic shifts depend on both 
V — I and V . At faint V magnitudes, this dependence of 
the shifts on both V and V — I can be explained by the 
existence of a magnitude limit in both the V and I filters. 
At brighter magnitudes, however, there remains a tendency 
for blue stars to be recovered redder than their true colours. 
This is likely an effect of crowding with the predominantly 
red cluster stars in the PC frame. 

Our WFPC2 frames are a mix of cluster and back- 
ground field stars, with no portion of the images contain- 
ing only one type of star. However, because the cluster and 
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Figure 10. Systematic and random photometric errors calcu- 
lated through artificial star tests for the example of NGC 1754. 
In (a), the arrows show the mean direction and distance over 
which stars originating in the bin at the tail of the arrow drifted 
in the colour-magnitude plane. In (6), the error bars show the la 
robust standard deviations of the recovered colours and magni- 
tudes of the artificial stars in each bin. For similar plots for the 
other clusters, see Olsen (1998). 



field stars are distributed differently, we can still statistically 
separate the two populations. We are further aided by the 
fact that the PC frames contain more cluster stars than field 
stars, vice versa in the WF frames. Our approach is to use 
the WF frames to remove field stars from the PC frame, 
and use the cleaned PC frame to remove cluster stars from 
the WF frames. After a few iterations of this procedure, we 
achieve satisfactorily clean cluster and field CMDs. 
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Our cleaning procedure accounts for the different 
crowding conditions and the different areas subtended by 
the cluster and field star distributions in the PC and WF 
frames. The following equation describes the number of con- 
taminating background stars expected in the frame contain- 
ing the stars of interest: 



N{V, V-I) 



APjV, V-I)J F{r)dA^ 
A=(F,y-7) jF{r)dAc 



D^{V,V-I) (1) 



where A^{V, V — I) is the crowding surface of the primary 

frame, A'^(V, V — I) is the crowding surface of the back- 
ground frame, -F(r) is the spatial distribution of the back- 
ground stars, dA^ and dA'^ are area elements, and D'^ is the 
distribution of stars in the CMD in the background frame. 
We assume the field stars are uniformly distributed, so that 
F{r) =constant and the integrals reduce to the areas of the 
frames. For the cluster stars, we choose a King (1966) model 
to represent F{r). The next section discusses these King 
model fits. 



4.3 Cluster profiles 

We determined the centers of the clusters through an al- 
gorithm developed by King and described in Mateo (1987). 
Briefly, the algorithm finds the center of symmetry of the 
clusters by extracting subimages on a grid of points in the 
image containing the cluster and correlating each subimagc 
with itself, rotated 180 degrees. The value of the correlation 
is stored at each grid point, resulting in a surface whose max- 
imum represents the center of symmetry. We measured the 
maximum by fitting this surface to a second order polyno- 
mial and computing the vertex. The algorithm worked well 
for the rich clusters, but was somewhat time-consuming and 
more difficult in the poorer clusters, particularly NGC 1898. 



5 CLEANING BACKGROUND STARS FROM 
THE CMDS 

With the cluster centers established, wo calculated the ra- 
dial stellar distributions of each cluster. We transformed the 
coordinates of all the stars to a reference frame centered on 
the cluster and scaled to the size of PC pixels. We set up 
thirty logarithmically spaced annuli around the center of 
each cluster out to the edges of the WF chips and calcu- 
lated the area of the portion of each annulus covered by the 
WFPC2 chips. We then summed up the number of stars 
within each annulus, divided the sum by the completeness 
appropriate to the stars' positions and magnitudes, and di- 
vided by the appropriate area. Because the magnitude limit 
is a function of position, we tried to restrict the radial dis- 
tribution to contain only stars with completeness ^ 50 per 
cent at all radii. In order to have enough stars to define the 
profiles of NGC 1835, NGC 2005, and NGC 2019, however, 
we were forced to include stars with completenesses as low 
as 20 per cent. 

We fit King (1966) models with variable amounts of 
background to the radial distributions to determine the rel- 
ative contribution of the cluster stars in each chip. First, 
we estimated the peak density from the profile close to the 
center and set the core radius, rc to the radius at which the 




Figure 11. By-eye fits of King (1966) model profiles (solid lines) 
to the completeness-corrected stellar surface density profiles of 
each cluster (open squares). The error bars are due to counting 
statistics and uncertainty in the completeness corrections. 



density reaches ~ 1/2 the peak radius. We then tried dif- 
ferent variations of r^, concentration, and background level 
and chose the best fit by eye. We note that for NGC 1835, 
the innermost three points do not give a good fit, probably 
because of incompleteness problems near the cluster center. 
Figs, lla-e show the observed background-subtracted pro- 
files for each cluster and the adopted fits. The parameters 
of these fits are included in Table 6. 



5.1 Removal of stars 

We began by removing a statistical sample of field stars 
chosen from the combined WF CMDs from the cluster- 
dominated PC CMD. Equation 1 becomes: 



iV^'^(F,V-7) 



A^^(y,y- J) A 



PC 



-D^^(V,y-7) 



AWF{V,V -I) A'^P 

N^'^iy, V - I) represents how the stars in the WF CMD 
would be distributed in the PC CMD, and are the stars 
we wished to subtract. However, because some areas of the 
CMDs are sparsely populated, not all of the bins in the 
PC CMD for which N^'^iy, 1/ - 7) is positive contain stars. 
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Figure 12. Color-magnitude diagrams of the LMC clusters after selection of stars of DoPHOT object type 1 with r > 9" from the 
cluster centers. For all of the clusters except NGC 1916, field stars have been statistically removed. The error bars are those calculated 
from the artificial star tests discussed in the text. 



Therefore, rather than simply subtracting the Af^'" (F, V—I) 
stars from the appropriate bins, we subtracted stars based 
on the error eUipse at each point. At each point (V, V — I), 
we gave the surrounding stars weights based on the error dis- 
tribution at {V, V — I), and chose the stars to be subtracted 
randomly from this weighted distribution. Stars outside the 
3-cr error ellipse were given zero weight and were not con- 
sidered. If no stars laid within the 3-a ellipse, we subtracted 
the nearest available neighbor. We then used the cleaned 
PC CMDs to subtract cluster stars from the WF CMDs. 
The relevant equation is: 



KPC(y^V-I) J FKir.g{r)dAPO 



x£>^^(y,l/-7) 

where JkingC?") is the King model fit from the completeness- 
corrected stellar density distribution and D^'^ {V, V — /) is 
the cleaned PC CMD. We used the same procedure de- 
scribed above to remove stars from the WF CMDs, pro- 
ducing CMDs with mostly cluster stars removed. 
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Because there are both cluster and field stars in both 
the PC and WF CMDs, we necessarily removed some clus- 
ter stars from the PC CMDs and some field stars from the 
WF CMDs in the first iteration of the subtraction procedure. 
Repeated iterations improved the separation, producing sta- 
tistically cleaner cluster and field star CMDs. We found that 
the procedure converges in ~5 iterations, as expected by the 
measured contamination levels and a simple calculation of 
the convergence rate. 



6 ANALYSIS OF CLUSTER CMDS 
6.1 Overview 

Figs. 12a-c show the cleaned cluster CMDs. Most field stars 
have been removed, although a few remain between the 
MSTO and the HB. Overall, however, the cleaning process 
makes it much easier to see the location of the MSTO. The 
CMDs all clearly show that the clusters are old, with MSTOs 
at V~23 and blue HB stars along with red ones. Compar- 
ing our CMD of NGC 1754 with the ground-based CMD 
of Jensen, Mould, & Reid (1988), we see the immense dif- 
ference that the resolution of HST mates in these crowded 
fields. While Jensen et al.'s (1988) effort was valiant, their 
misidentification of the field main sequence and the HB of 
NGC 1754 with the cluster main sequence led them to sug- 
gest an age of only 0.8 Gyr for the cluster. 

It is tempting to suggest that some of the blue stars just 
above the MSTO in NGC 1898 and NGC 2019 may be true 
blue stragglers. However, considering the large numbers of 
young MS field stars in this region and the uncertainty in 
the statistical cleaning procedure, we malce no such claim 
here. With the exception of NGC 2005 and possibly NGC 
1916 (Fig. 7d), all of the clusters have HB stars both to the 
blue and to the red of the instability strip. Both the NGC 
1835 and NGC 2019 CMDs show a few AGB stars, as do 
possibly the other clusters. 

NGC 1898 is unique in containing a large number of 
stars to the red of the primary fiducial secjucnce. As we were 
puzzled by these stars, we selected them on the CMD and 
examined their locations in the image. Figs. 13a&b show 
the selected stars and their locations on the image. The ma- 
jority of the stars lie in an area of the sky that contains 
relatively few stars-this plus the anomalously red colours of 
these stars indicates that the area is a patch of heavy ob- 
scuration. Amazingly, the rest of the cluster suffers only a 
modest of amount of reddening, comparable to that of the 
other clusters. 
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Figure 13. (a) The peculiarly red stars in the colour-magnitude 
diagram of NGC 1898 were selected within the bounding polygon 
shown and their positions displayed on the image in (6). The 
selected stars are heavily concentrated towards a region of the 
image with noticeably lower star density, and are therefore very 
likely highly reddened. 



6.2 Comparison of LMC and Milky Way globular 

cluster systems 

We used our knowledge of the abundances of the LMC 
clusters in conjunction with available methods to measure 
the ages of the clusters. Due to the strong dependence on 
the uncertain photometric zero points, we avoided measur- 
ing the ages from a direct comparison of the CMDs to 
model isochrones. Instead, we relied on differential tech- 
niques which axe free of zero point dependence. We used the 
approach of VBS to compare our CMDs to fiducial sequences 



of Milky Way globular clusters of similar abundance, mea- 
suring relative ages with an accuracy of < ±1 Gyr from the 
difference in position of the RGBs. We also measured AI/Jjb, 
which is a calibrated function of age. While less accurate 
than the "horizontal" VBS method, AT/hb has the advan- 
tage of being fundamentally better understood. Finally, as 
illustrated by LDZ, comparison of cluster HB morphologies 
with HB evolutionary tracks is a potentially very powerful 
way to establish the chronology of globular cluster formar 
tion. While it is controversial, we used this technique to at- 
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Figure 14. The colour-magnitude diagrams of the LMC clusters shifted to match Milky Way globular cluster fiducial sequences in 
turnoff colour and HB magnitude. The dashed line to the blue of the fiducial RGB shows the location it would have if the fiducial were 
2 Gyr older, while the line to the red shows the RGB of a cluster 2 Gyr younger than the fiducial. The locations of these lines were 
calculated from the models of VandenBerg (1997). 



tempt to build a consistent picture of the LMC cluster ages. 
Where inconsistencies arise, we can then interpret them to 
be due either to errors in the observations or a failure in 
some portion of the theoretical framework. 

Each of the above age-dating techniques depends on the 
assumed metallicity, the VBS technique through the need 

to compare clusters of similar metallicity, AV^^g through 
the metallicity dependence of both the magnitude of the 
HB and the magnitude of the MSTO, and the LDZ tech- 
nique through metallicity being the "first parameter" af- 
fecting HB morphology. While we have available directly 
measured spectroscopic abundances of the clusters (091), 
these abundances are based on difficult observations of only 
1-2 giant stars per cluster. We therefore measured the abun- 
dance of each cluster from the CMD using the formalism of 
Sarajedini (1994, hereafter S94), described in detail below. 
Because these CMD-basod abundances arc rooted in well- 
understood photometric data of a large number of stars, we 
regard them as more reliable than the currently available 
spectroscopic abundances. 



6.2.1 Abundances and reddenings from the CMDs 

Wc followed S94 to derive abundances from the height of the 
RGB above the level of the HB; this method also yields the 
reddening from the colour of the RGB at the level of the HB. 
As suggested in S94, we fit the portion of the RGB brighter 
than the HB with a 2nd order polynomial, which let us solve 
for E{B — V) and [Fe/H] analytically. We determined the 
level of the HB by taking the median of the points on the flat 
portion of the HB. For NGC 2005, which has mostly blue 
HB stars, this involves some guesswork, so its parameters 
are likely more uncertain than those of the other clusters. 

To calculate the uncertainties, we performed the poly- 
nomial fits to the RGB on a set of 100 Monte Carlo re- 
alizations of each CMD. These simulated CMDs were con- 
structed by choosing the appropriate number of stars from 
the luminosity functions of Vandcnbcrg (1997) isochrones 
with abundances and reddenings approximately matching 
the observed CMDs. We simulated the observational errors 
by applying the V — I and V shifts calculated from our 
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Figure 14. cont. 



artificial star database. After introducing a 0.1 magnitude 
Gaussian-distributed error in the V magnitude of the HB, 
which was not simulated, we then calculated E{B — V) and 
[Fe/H] for the simulated CMDs and compared the values 
to the input values. We found typical errors of 0.15 dex in 
[Fe/H] and 0.01 magnitudes in E{B — V) with systematic 
biases of up to 0.1 dex in [Fe/H] and 0.025 magnitudes in 



E{B — V). These biases are probably due to small scale 
changes in the photometry introduced by crowding, as illus- 
trated in the example of Fig. 10. Table 6 contains the cal- 
culated E(B — V) and [Fe/H] values for the clusters, after 
correction for the biases found in the Monte Carlo experi- 
ments. 
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Figure 15. Calibration of age difference as a function of colour 
offset of the RGB compared to the calibration of Harris et al. 
(1997) for the example of [Fe/H]= -2.14. The caUbration differ- 
ences are small compared to the precision of the age measure- 
ments. 

6.2.2 Relative ages through comparison with Milky Way 

clusters 

We compared our LMC cluster CMDs to fiducial sequences 
of M3, M5, M13, M92 (Johnson & Bolte 1997), and M55 
(Manduslicv et al. 1996 and personal communication), which 
span the range of metallicity —1.4 <^ [Fe/H] <^ —2.2. For 
consistency, we adopt the Zinn & West (1984) abundances 
for all of the Milky Way clusters. We used the technique 
described by VBS to establish the comparison between the 
LMC clusters and the Milky Way fiducials. With this tech- 
nique, the cluster CMD and a fiducial of similar metallicity 
are registered in the C-M plane through the colour of the 
main sequence turnofF (MSTO) and the magnitude of the 
point 0.05 magnitudes redder than the MSTO. The diflfer- 
cnce in the lengths of the subgiant branches can then be 
easily measured and converted to a relative age tiirough a 
calibration from model isochrones. 

As an alternative to the process described above, which 
is fairly automatic and unbiased, we also attempted to de- 
rive comparisons by eye, involving both the disadvantage of 
possible subconscious bias and the advantage of judgement 
based on knowledge. Each of the six LMC cluster CMDs 
was compared to the fiducials of M3 and MIS and displace- 
ments inV — I and in V were recorded for the best eye fits. 
(This was done by an author who had not, up to that time, 
been involved with the initial fitting process). Interestingly, 
the results were quite similar to those from the automatic 
comparison. The average difference in the derived V — I was 
0.01 magnitudes, while that in V was 0.11 mag. These by- 
eye measurements confirm that the automatic method had 
not introduced any unreahstic data into the problem. 



Table 5. Calibration of Relative Ages 
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6.2.3 Registration of the CMDs 

We measured the colour of the MSTO, (V — /)msto, of 
our clusters and the fiducials by fitting a second order poly- 
nomial to a region around the turnoff and calculating the 
bluest point of the fit. We iterated the fit until the V — I 

colour of the bluest point of the polynomial fit changed by 
less than 5x10"* magnitudes, each time rejecting outlying 
points. To measure the V magnitude of the point 0.05 mag- 
nitudes redder than the MSTO (1^+0.05) wc used a similar 
fitting procedure. We fit a straight line to points near V+0.05, 
rejecting outliers and iterating the fit until V+0.05 changes 
by less than 0.01 magnitudes. However, we found that af- 
ter using V+0.05 to register the LMC cluster CMDs with the 
Milky Way fiducials, the magnitudes of the HBs occasionally 
disagreed by > 0.1 magnitudes, even for clusters of the same 
metallicity. In these cases, we consider our eye estimates of 
the level of the HBs more reliable than the measured V+o.os 
points, as these V+0.06 measurements are likely affected by 
incompleteness in the photometry. 

Figs. 14a-i show the comparison of our LMC clusters 
with Milky Way fiducials after registration. The dashed lines 
indicate the differences in position of the RGB that would 
be expected for a 2 Gyr age difference. We only show those 
comparisons for which the abundances of the fiducials most 
closely match either the spectroscopic or CMD-based abun- 
dances. For NGC 1754, NGC 1835, and NGC 1898, the dif- 
ference between the comparisons based on the spectroscopic 
abundances and those based on the CMD abundances are 
small. For NGC 2005 and NGC 2019, the CMD-based abun- 
dances are ~0.6 dex higher than the spectroscopic ones. 
The high CMD-based abundance of NGC 2005 especially 
comes as a surprise because the cluster has a very blue HB 
morphology. While the discrepancy with the spectroscopic 
abundances is a cause for some concern, wc note that the 
spectroscopic abundances of NGC 2005 and NGC 2019 are 
based solely on measurements of a single star in each cluster, 
which could be in error. In addition, the CMD-based abun- 
dances produce better matches of the Galactic fiducials to 
the overall CMDs, reinforcing our greater confidence in these 
abundances. 

6.2.4 Relative ages and age errors 

To measure the difference in position of the RGBs in Figs. 
14a-i, wc fit the piece of the fiducial extending from —5 < 
AV < -2 to the CMD points in the same AT/ range. To 
map the differences in RGB position to age differences be- 
tween the Galactic and LMC clusters, we used a calibration 
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based on the new isochrones of VandenBerg (1997). These 
isochrones are calculated for more than a dozen metallicities 
in the range -2.3 < [Fe/H] < -0.3 for ages of 8-18 Gyr in 
2-Gyr steps and fit observations of M92 extremely well. Af- 
ter registering the six isochrones of a given mctallicity in the 
A{V — I), AV plane, we fit a piece of the 14 Gyr isochrone 
spanning the range —5 < AV < —2 to similar pieces of 
the other five isochrones. To establish the calibration over 
a continuum of age differences, we fit a 2nd order polyno- 
mial to the offsets needed to shift this portion of the RGB of 
the 14 Gyr isochrone to the RGBs of the other isochrones. 
The choice of the 14 Gyr isochrone as the fiducial for this 
cahbration is inconsequential, as the change of — / coloxir 
of the RGB is nearly linear with age for ages between 12 
and 18 Gyr, the likely age range of old globular clusters. 
Table 5 lists the polynomial coefficients of the fits to the 
equation AAge = a + bA{V -I) + cA{V - 7)^ for the range 
of metallicities considered in this paper. Fig. 15 shows our 
calibration for [Fe/H]=— 2.14 compared to that of Harris et 
al. (1997). Over a range of age differences of ±1 Gyr, our 
calibration differs by ~15 per cent from that of Harris et al. 
(1997). This difi^erence may be attributed to the fact that 
wo used a large portion of the RGB to measure the age dif- 
ferences while Harris et al. (1997) apparently used only the 
difference in RGB color at AV = 2.5. Because we can more 
accurately measure the differences in RGB colors if we use 
many points on the RGB, we chose to use our calibration 
over that of Harris et al. (1997). However, as the difference 
in the calibration will only produce differences in the rel- 
ative age measurements of <0.2 Gyr, the choice has little 
consequence for our results. 

Table 6 contains the ages we derive, in Gyr, for the LMC 
clusters relative to the appropriate comparison clusters. On 
average, the LMC clusters have the same ages as the Milky 
Way clusters to within 1.0 ± 1.3 Gyr, reinforcing our ear- 
lier statement that the LMC old globular clusters arc very 
similar to those in the Milky Way. In order to examine this 
global statement on a cluster-to-cluster basis, however, we 
need to establish the accuracy of the individual relative age 
measurements. We did this by performing the relative age 
measurement on the sets of 100 Monte Carlo realizations de- 
scribed previously. For each simulated CMD, we measured 
{V — /)msto and V^+o.o5 as for the observed CMD and regis- 
tered the simulated CMD and the error-free input isochrone 
in the A{V — I), AV plane. The simulated measurements of 
V+0.05 confirm that low completeness affects this measure- 
ment, justifying our use of the HB V magnitude to perform 
the vertical registration. Since we do not simulate the HB 
stars, we assumed that we can correctly perform the vertical 
registration to within 0.1 V magnitudes. Thus, we vertically 
registered the simulated CMD according to the exact V+o.os 
value calculated from the isochrone, but applied a random 
offset selected from a Gaussian with a 1-cr dispersion of 0.1 
magnitudes. We then measured the age difference between 
the simulated CMD and the input isochrone as we did for the 
real clusters. Figs. IGa^e show the resulting distributions of 
age differences and the cumulative distributions of absolute 
age differences around the mean. The distributions of age 
differences are narrow, with 70 per cent of the distributions 
enclosed within 0.7-0.8 Gyr of the mean. However, three of 
the clusters, NGC 1754, NGC 2005, and NGC 2019, show 
systematic biases towards recovering ages ~1 Gyr younger 



than the age of the input isochrone, while NGC 1835 has a 
bias towards recovering ages ~0.5 Gyr younger than the in- 
put. These biases are likely due to the crowding errors in the 
photometry illustrated by the example of Fig. 10. Although 
ignoring the biases would not significantly affect our results, 
we have corrected for them in the relative age measurements 
reported in Table 6. 

6.2.5 Ages from AV'i^ 

To measure AV^^g we used the values of Vhb measured in 
section 6.2.1 and the magnitude of the MSTO (Vmsto) im- 
plied by the polynomial fit to the points near the turnoff. 
Because the MSTO region is nearly vertical in the CMD and 
the measurement of Vhb is subject to some interpretation, 
the error in AV'Jj^ may be larger than reported here. We 
estimate the error in Vhb to be ±0.1 magnitudes, while for 
the error in Vmsto we adopt the error calculated from the 
covariance matrix of the polynomial fit. 

To convert AV^hb to age, we adopted the calibration of 
Chaboycr, Dcmarque, & Sarajedini (1996), which is based on 
a recent set of Yale isochrones (Chaboyer & Kim 1995). The 
calibration depends on an assumed Mv(RR) — [Fe/H] rela- 
tion; their preferred relation is Mv(RR) = 0.20[Fc/H] +0.98. 
Recent work based on HIPPARCOS parallaxes and proper 
motions of field RR Lyrae (Tsujimoto, Miyamoto, & Yoshii 
1998) finds a similar slope and zero point, Mv(RR) = 
0.20[Fe/H] -I- 0.91, albeit with large error. The indepen- 
dent analysis of Fernlcy ct al. (1998) gives Mv{KR) = 
0.18[Fe/H] -I- 1.05, which agrees, within the error, with the 
Tsujimoto et al. (1998) result. As the Chaboyer et al. (1996) 
preferred zero point of the My (RR) — [Fe/H] relation is the 
exact average of the two available HIPPARCOS-hased zero 
points, we adopt their preferred relation. An advantage of 
using this relation is that we can easily compare the ages 
we derive from AI/Jb with those of the Galactic globu- 
lar clusters. As has been noted, however, this relation dis- 
agrees strongly with the distance to the LMC implied by the 
first HIPPARCOS Cophcid work (Feast & Catchpole 1997), 
but not with the recent HIPPARCOS Cepheid analysis of 
Madore & Preedman (1998). The AV^b ages we measure 
arc listed in Table 6, for both the 091 and CMD-based 
abundances. 



6.2.6 Analysis of horizontal branch morphology 

An interesting question is whether the LMC clusters show a 

variation of HB morphology with mctallicity similar to the 
Milky Way clusters. We measured the HB morphology of 
our clusters through the connuouly used (B-R)/(B-|-V-|-R) 
index, where B is the number of stars to the blue of the insta- 
bility strip, R the number of stars to the red of the instability 
strip, and V the number RR Lyrae variables. Because we do 
not have sufficient time resolution in our images to identify 
RR Lyrae variables from their light curves, we rely on the 
approximate location of the instability strip in the CMD 
plane to determine the boundaries of the B, V, and R zones. 
Our chosen boundaries are V — 1—0.23 for the blue edge of 
the instability strip and V — 7=0.57 for the red edge. Ta- 
ble 6 contains our measured values of (B-R)/(B-|-V-|-R) for 
each cluster, corrected for incompleteness, along with uncer- 
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Figure 16. Plots of the distribution of errors in the relative age measurement and cumulative distribution of errors around the 
calculated from the Monte Carlo experiments described in the text. 70 per cent of the cumulative distributions are enclosed within 
Gyr. However, each cluster shows a bias of 0-1 Gyr towards measuring younger relative ages. 
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taintios calculated from Poissori counting statistics and the 
uncertainties in the completeness corrections. 

In Figs. 17a&b we plot [Fe/H] vs. (B-R)/(B+V+R) for 
the LMC clusters studied here, for the outer LMC clusters 
using data from Walker (1992b), and for selected Milky Way 
clusters from Table 1 of LDZ. We show separate plots for the 
091 spectroscopic abundances and the CMD-based abun- 
dances. Depending on which set of abundance measurements 
is adopted, we arrive at slightly different conclusions. If we 
choose to adopt the spectroscopic abundances, we conclude 
that, based on the HB morphologies, our LMC clusters are 
similar to the "younger halo" of the Milky Way. On the other 
hand, if we adopt the CMD-based abundances, we find that 
the LMC clusters are similar to the oldest Milky Way glob- 
ulars, with NGC 1835 perhaps being ~2 Gyr younger. As 
we are more confident in the CMD-based abundances, the 
HB models suggest that the LMC is as old as the Milky 
Way. Independent of the abundances used, all of the clus- 
ters, with the possible exception of NGC 1835, fall on the 
same HB evolutionary sequence. The similarity in age of the 
clusters implied by the HB sequences agrees with the narrow 
age spread implied by both the comparison with Milky Way 
fiducials and the measurements of AF^b. It is also inter- 
esting to note that none of the LMC clusters studied here 
bears similarity to Ruprecht 106 or Pal 12, which have been 



suggested to be captures from the Magellanic Clouds (Lin 
& Richer 1992; cf. Rup 106 CMD of Buonanno et al. 1993). 



6.2.7 Summary of reddenings and distances 

We estimated the reddenings of the clusters in two ways. 
First, the horizontal shifts needed to register our CMDs with 
the the Milky Way fiducials yields the reddening difference 
between the LMC clusters and the Milky Way comparison 
clusters. As the reddenings of the Milky Way clusters is well 
known, we can calculate the reddenings of the LMC clusters. 
Second, as previously discussed, we use the S94 method to 
measure the reddening from [V — 7)g, the colour of the giant 
branch at the level of the HB. Table 6 contains our measure- 
ments of E{B — y) as well as the estimates of Suntzeff et al. 
(1992). For the Milky Way comparison method, we report 
separate values for the two possible abundance systems. In 
all calculations, we have used E{B — V) = E{V — 7)/1.3 
(Dean, Warren, & Cousins 1978). As the dispersion in the 
measurements is small, we adopt for the E{B — V) of each 
cluster the average of our measurements, also listed in Table 
6. 

Our CMDs are not sufficiently deep to accurately mea- 
sure the distances to the clusters directly from a fit to the 
unevolved main sequence. However, we can use the vertical 
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Table 6. Derived Parameters of LMC Globular Clusters 
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[Fe/Hl 


091 


-1.54 ± 0.2 


-1.79 ± 0.2 
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E(B-V) 
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0.08 ±0.02 


0.13 ± 0.02 
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09 1 abundance 
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0.12± 0.02 
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15.6 ± 2.2 


16.2 ± 2.8 


13.5 ± 2.2 
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18.60 ± 0.16 


18.65 ± 0.16 


18.45 ± 0.16 


18.40 ± 0.16 
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18.60 ± 0.12 


18.40 ± 0.12 
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shifts measured in Section 6.2.3 to measure the distances 
to the LMC clusters based on adopted distances to the 
Milky Way comparison clusters. In addition, by adopting 
an Mv(RR) — [Fe/H] relation, we can measure the distances 
to the LMC clusters from Vhb- By demanding that the dis- 
tances measured with our data should be consistent with 
those implied by observations of LMC field RR Lyrae, we 
can use the distance measurements to check the internal con- 
sistency of our measurements of the reddenings and abun- 
dances, in addition to establishing the relative distances be- 
tween the clusters. For our calculations, we adopt distances 
to M3, M5, and M55 that arc consistent with the LMC 
modulus of 18.5 used in Walker's (1992c) analysis of the 
zero point of the My (RR) — [Fe/H] relation; we use the pre- 
ferred Mv{RR) - [Fe/H] relation of Chaboyer et al. (1996), 
Mv(RR) = 0.20[Fe/H] ± 0.98, which is consistent with the 
SN1987A distance to the LMC of 18.37 (Gould 1995). While 
other analyses of the SN1987A ring find larger LMC dis- 
tances (Panagia et al. 1991, Sonneborn et al. 1997), the ex- 



act distance is irrelevant for the consideration of the internal 
consistency of our measurements. 

Fig. 18 shows our distance modulus measurements for 
both the Milky Way comparison and Vhb methods and for 
both the 091 and CMD-based abundance sources. We com- 
pare our measurements with the distance modulus gradient 
across the LMC of Caldwell & Coulson (1986). In the top 
panels, (a) and (b), we have used the 091 abundances while 
in the lower panels, (c) and (d), we have used the CMD- 
based abundances. The panels on the left, (a) and (c), show 
the results derived from the comparison with Milky Way 
clusters while the panels on the right, (b) and (d), show 
the results of adopting Mv(RR) = 0.20[Fe/H] ± 0.98 in con- 
junction with Vhb. We have subtracted an LMC distance 
modulus of 18.5 from the left panels and 18.37 from the 
right panels. While the errors are large, the plots suggest 
that the CMD-based abundances produce a greater consis- 
tency between the cluster distances and the LMC modulus. 
All of the distances, with the possible exception of NGC 
1754, appear consistent with lying in the plane of the LMC. 



© 0000 RAS, MNRAS 000, 000-000 



0.0 |—r 



-2.5 1 I , 
-1.0 



[Fe/H] from RGB 



-0.5 



0.0 

(B-R)/(B+V+R) 



0.5 



(a) 




0.0 



-0.5 



-1.0 



Old Globular Clusters in the LMC 21 

Spectroscopic [Fe/H] 



X 



-1.5 



-2.0 



-2.5 



1.0 



-1.0 



-0.5 



0.0 

(B-R)/(B+V+R) 



0.5 



(b) 




1.0 



Figure 17. Compajrison of the horizontal branch morphologies of the LMC clusters of this study (filled circles), LMC clusters from 
Walker (1992b) (open circles), and Milky Way comparison clusters used in this study (crosses) using the the CMD-based abundances 
(a) and the spectroscopic abundances of Olszewski et al. (1991) (6). The solid lines are HB evolutionary tracks from Lee et al. (1994) 
showing age differences of ±2 Gyr. 



7 SUMMARY AND CONCLUSIONS 

The main results of our analysis are that the six old LMC 
globular clusters NGC 1754, NGC 1835, NGC 1898, NGC 
1916, NGC 2005, and NGC 2019 arc clearly very similar in 
age, abundance, and HB morphology to the globular clus- 
ters of the Milky Way halo. Excluding NGC 1916 from the 
bulk of the analysis because of the difficulty of treating its 
differential reddening, we have explored the similarity in age 
through the use of three age-dating techniques: comparison 
of our CMDs with Milky Way globular cluster fiducial se- 
quences through the "horizontal method" (VBS), measure- 
ment of Al^HB ("vertical method"), and comparison of our 
HB data with HB evolutionary models from LDZ. While 
showing the clear similarity in age with Milky Way globular 
clusters, each of these techniques additionally suggest that 
we detect no internal age spread in this set of LMC globular 
clusters. We have measured the abundances of the clusters 
from the slope of the RGB (S94). While the mean abun- 
dance that we derive is higher than that measured by 091, 
it is not grossly different. We have argued, however, that the 
CMD-based abundances should be considered more reliable. 

Throughout our analysis, we have tried to be as careful 
as possible in identifying and measuring errors. In particular, 
we have extensively used Monte Carlo simulations to model 
both systematic and random errors in our measurements. 
Although we were unable to independently check the Holtz- 
man (1995b) zero points to a precision of <0.1 magnitudes 
through our ground-based photometry, our meEisurements 
of the abundances and relative ages are differential, and so 
should not be affected by errors in the zero points. The red- 
denings and distance moduli that we derive, however, are 
subject to unkown zero point errors in the photometry. 

Adopting the CMD-based abundances, the possibility 



that the old LMC globular clusters are younger than the 
oldest Milky Way globular clusters by 2-3 Gyr (Da Costa 
1993) appears not to be the case for the clusters studied 
here. This result is clear from the relative ages derived from 
the comparison of the clusters with Milky Way fiducials, the 
resuhs of which are shown in Table 6. NGC 1754, NGC 1898, 
NGC 2005, and NGC 2019, which we find have abundances 
most closely matching that of the classical old Galactic glob- 
ular cluster M5, are all of similar age or older than M5. NGC 
1835, which we find has an abundance similar to M3, is ~2 
Gyr older. As LDZ consider M3 to be ~2 Gyr younger than 
the oldest halo clusters, NGC 1835 appears to have an age 
similar to the oldest Milky Way clusters. 

Interpreting our data through the HB models of LDZ, 
we paint a similar picture. As shown in Fig. 17a, the clusters 
NGC 1754, NGC 1898, NGC 2005, and NGC 2019 all fall on 
the HB evolutionary track crossing through the oldest Milky 
Way halo clusters. Comparing our data with the HB data of 
Walker (1992b) for the outer LMC old globular clusters, we 
find that the HB models suggest that the earliest episode of 
cluster formation in the LMC spanned 2-3 Gyr. However, we 
find some problems with the ages implied by the HB models. 
The age of NGC 1835 suggested by the HB models is incon- 
sistent with its age relative to M3. In addition, the age of 
NGC 1754 relative to M5 is slightly older than suggested by 
the HB models. These inconsistencies may indicate a diffi- 
culty with the interpretation that mctallicity and age arc the 
only parameters affecting HB morphology, as has been sug- 
gested by many authors (e.g. Buonanno et al. 1997, Catelan 
et al. 1997, Sweigart 1997). Given the currently large errors 
in the abundance measurements, however, we do not claim 
on the basis of our data that this interpretation needs to be 
modified. 

The ages we derive from AT/hb are marginally consis- 
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Figure 18. Comparison of the distance moduli of the LMC clus- 
ters implied by two different methods using two sources of abun- 
dances with the model of the tilt of the plane of the LMC of 
Caldwell & Coulson (1986). (o) and (b) incorporate the spectro- 
scopic abundances of Olszewski et al. (1991) while (c) and (d) 
use the CMD-based abundances, (a) and (c) show the distances 
implied by the comparison to Milky Way clusters while (b) and 
(d) show the moduli implied by the horizontal branches, assuming 
Mv(RR) = 0.20[Fc/H]-F0.98 (Chaboycr ct al. (1996). A modulus 
of 18.5 has been subtracted from (a) and (c) and 18.37 from (6) 
and (d) to show the consistency with the distances implied by 
LMC RR Lyrae observations and the adopted zero points of the 
Mv'(RR) - [Fe/H] relation. 



tent with the picture that the LMC clusters studied here 
are as old as the oldest Milky Way clusters. Our measured 
values of AI/^b imply an average age of 15.3 ±1.5 Gyr, 
which is lower but within ~ 1.5a of the average age of the 
old halo quoted by Chaboyer et al. (1996), 17.8 ± 0.4 Gyr. 

Because of the metallicity dependence of the age in- 
dicators discussed in Section 6.2, our conclusions based on 
the cluster ages depend critically on the CMD-based abun- 
dances being correct. If we adopt the 091 spectroscopic 
abundances, for instance, our conclusions change consider- 
ably. For example. Fig. 17b shows the comparison of the 
LDZ HB evolutionary tracks with our cluster data using the 
091 abundances. The HB tracks imply, in this case, that 
the LMC clusters are indeed ~2 Gyr younger than the old 
Milky Way halo. However, we are then unable to build a 
consistent picture with the relative ages derived from the 
comparison to Milky Way clusters. Coupled with the better 
internal consistency weakly implied by the distance moduli 
(Figs. 18a-d), we continue to prefer the CMD-based abun- 
dances for these LMC clusters. Nevertheless, it would clearly 
be extremely valuable to have high-resolution spectroscopic 
abundances of several stars in each cluster available. 

In contrast to the globular clusters of the Milky Way 
halo, the old LMC clusters are not as clearly part of a dis- 
tinct halo (091; Schommer et al. 1992). Our results imply 



that at the time the Milky Way formed its first globular 
clusters, the LMC may have already collapsed to a disk and 
started forming clusters. Why there appears to be no clear 
halo component in the LMC (Olszewski, Suntzeff, & Mateo 
1996) remains an open question. However, by establishing 
that the oldest clusters in the LMC are truly as old as the 
Milky Way Galaxy, we have taken am important step towards 
understanding the formation of the LMC. 



ACKNOWLEDGEMENTS 

We arc indebted to Michael Bolte and Jennifer Johnson for 
providing us with the unpublished {V — I),V fiducial se- 
quences of M3, M5, MIS, and M92 and to Georgi Mandushev 
for the unpublished sequence of M55. We are grateful to Don 
VandenBerg for his invaluable help with the cluster ages and 
to Carlton Pryor for his help with the King model fits. We 
thank Eric Deutsch for his modifications to DoPHOT which 
made it easier to use. The anonymous referee is thanked 
for helpful comments which improved the manuscript, as is 
George Wallerstein for his reading and comments. This work 
was supported by STScI grant GO05916. 



REFERENCES 

Brocato E., Castellani V., Ferraro F. R., Piersimoni A. M., Testa 

v., 1996, MNRAS, 282, 614 
Buonanno R., Corsi C. E., Pecci F. F., Richer H. B., Fahlman G. 

G., 1993, AJ, 105, 184 
Buonanno R., Corsi C, BcUazzini M., Ferraro F. R., Pccci F. F., 

1997, AJ, 113, 706 
Caldwell J. A. R., Coulson I. M., 1986, MNRAS, 218, 223 
Catclan M., Swcigart A. V., Borissova .)., Spassova N., 1997, in A 
Half Century of Stellar Pulsation Interpretations; A Tribute 
to Arthur N. Gox, eds. P. A. Bradley & J. A. Guzik, ASP 
Conference Series, vol. 135 
Chaboyer B., Kim Y.-C, 1995, ApJ, 454, 767 
Chaboycr B., Dcmarquc P., Sarajcdini A., 1996, ApJ, 459, 558 
Da Costa G. S., 1993, in The Globular Cluster-Galaxy Connec- 
tion, cds. G. H. Smith & J. P. Brodie, ASP Conference Series, 
vol. 48, p. 363 

Dean J. F., Warren P. R., Cousins A. W. J., 1978, MNRAS, 183, 
569 

Feast M. W., Catchpole R. M., 1997, MNRAS, 286, LI 

Fernley J., et al., 1998, A&A, 330, 515 

Gould A., 1995, ApJ, 452, 189 

Harris W. E. ct al., 1997, AJ, 114, 1030 

Hodge P. W., 1960, ApJ, 131, 351 

Holtzman J., ct al., 1995a, PASP, 107, 156 

Holtzman J., ct al., 1995b, PASP, 107, 1065 

Jensen J., Mould J., Rcid N., 1988, ApJS, 67, 77 

Johnson J. A., Bolte M., 1997, BAAS, 189, 7118 

King I. R., 1966, AJ, 71, 64 

Landolt A. U., 1992, A J, 104, 340 

Lee Y.-W., 1992, PASP, 104, 798 

Lcc Y.-W., Dcmarquc P., Zinn R., 1994, ApJ, 423, 248 (LDZ) 
Lin D. N. C, Richer H. B., 1992, ApJ, 388, L57 
Olscn K. A. G., 1998, Ph.D Thesis, University of Washington 
Olszewski E. W., Schommer R. A., Suntzeff N. B., Harris H. C, 

1991, AJ, 101, 515 (091) 
Olszewski E. W., Suntzeff N. B., Mateo M., 1996, ARA&A, 34, 

511 

Madore B. F., Freedman W. L., 1998, ApJ, 492, 110 



© 0000 RAS, MNRAS 000, 000-000 



Old Globular Clusters in the LMC 23 

Mandushev G. I., Fahlman G., Richer H. B., Thompson I. B., 

1996, AJ, 112, 1536 
Mateo M. L., 1987, Ph.D Thesis, University of Washington 
Mighell K. J., Rich R. M., Shara M., Fall S. M., 1996, AJ, 111, 

2314 

Panagia N., Gilmozzi R., Macchetto F., Adorf H.-M., Kirshner R. 

P., 1991, ApJ, 380, L23 
Reid N., Freedman W., 1994, MNRAS, 267, 821 
Sarajcdini A., 1994, AJ, 107, 618 (S94) 

Sarajedini A., Chaboycr B., Dcmarquc P., 1997, PASP, 109, 1321 
Schechter P. L., Mateo M., Saha A., 1993, PASP, 105, 1342 
Schommer R. A., Suntzeff N. B., Olszewski E. W., Harris H. C., 

1992, AJ, 103, 447 
Shapley H., 1930, Star Clusters, (New York: McGraw-Hill) 
Sonneborn G., Fransson C., Lundqvist P., Cassatella A., Gilmozzi 

R., Kirshner R. P., Panagia N., Wamsteker W., 1997, ApJ, 

477, 848 

Stetson P. B., VandcnBcrg D. A., Boltc M., 1996, PASP, 108, 560 
Suntzeff N. B., Schommer R. A., Olszewski E. W., Walker A. R., 

1992, AJ, 104, 1743 
Sweigart A. V., 1997, ApJ, 474, L23 

Testa v., Ferraro F. R., Brocato E., Castellani V., 1995, MNRAS, 
275, 454 

Tsujimoto T., Miyamoto M., Yoshii Y., 1998, ApJ, 492, 79 
VandenBerg D. A., Bolte M., Stetson P. B., 1990, AJ, 100, 445 

(VBS) 

VandenBerg D. A., 1997, personal communication 
Walker A. R., 1990, AJ, 100, 1532 
Walker A. R., 1992a, ApJ, 103, 1166 
Walker A. R., 1992b, AJ, 104, 1395 
Walker A. R., 1992c, ApJ, 390, L81 

Westerlund B. E., 1997, The Magellanic Clouds, Cambridge as- 
trophysics series, (Cambridge: Cambridge University Press) 

Whitmore B., Heyer I., 1997, WFPC2 Instrument Science Report 
97-08 

Zinn R., West M. J., 1984, ApJS, 55, 45 



© 0000 RAS, MNRAS 000, 000-000 



